LOCATION METHOD FOR MOBILE NETWORKS 

Field of the invention 
[01] The present invention relates generally to location techniques. More 

5 specifically, the present invention relates to determination of the geographical 

location of a mobile (i.e. a mobile terminal) within a mobile networlc. 



Baclcground of the invention 
[02] There are two major reasons that have given motivation and fueled the 

10 development of location determining techniques in mobile networks. First, dif- 
ferent authorities set requirements for the location determination of mobile 
terminals. It is highly desirable that certain authorities, such as emergency call 
centers, can locate the calling party as accurately as possible. In many coun- 
tries legislation sets requirements for such location methods. For example, in 

15 the USA mobile location to an accuracy of 50 meters for 67 percent of calls 
and 150 meters for 95 percent of calls will be mandatory in the near future for 
handset-based solutions. Second, many of the future services provided in 
mobile networks will be such that they require information about the current 
geographical location of the mobile terminal. 

20 [03] The greater the accuracy of a location method, the better it can serve 
the application utilizing the location information. This applies especially to 
densely built urban areas where long-range visibility is not possible. The accu- 
racy of the location methods is dependent on many different factors, such as 
radio propagation effects. This is a major reason why dense urban areas form 

25 a difficult and challenging environment in terms of location accuracy, with se- 
vere multipath characteristics and signal propagating with less attenuation 
along street canyons than through buildings. 

[04] When a mobile is located, successive measurement results are obtained 
from the network. The measurement results are typically not received as co- 
30 ordinates but as parameters which define a certain area within which the mo- 
bile is located with a certain probability. As the network normally provides 



several parameters fore single measurement, the parameters relating to a 
single measurement are in this context denoted as a parameter set. On the 
basis of the parameter set, geometrical limitations can be determined which 
define the geographical boundaries of one or more areas where the mobile is 
5 most likely located. As discussed below, the determination of the location in 
current cellular networks is widely based on a parameter set including the cell 
ID and the Timing Advance value, which indicate how far from a certain base 
station the mobile most probably is located. 

[05] The cunrent location methods employ the parameter sets obtained by 
10 first defining the said one or more areas and then calculating an estimate for 
the location of the mobile. History data, i.e. previous estimates, can then be 
used to make the current estimate more accurate. One known method, which 
O is based on a least squares scheme, determines a path for which the sum of 

the squares of the deviations from the estimated location points obtains a 
15 minimum value. 
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[061 A major drawback of the current methods is that they do not preserve 

the original data revealed by the parameters obtained from the mobile net- 

p work as to the geographical distribution of the location of the mobile. Instead, 

with current methods the estimation process moves in an early phase to a 

^ 20 point-based approach where the information concerning the geographical dis- 
cs 

tribution is reduced to estimated location points and the final estimate is calcu- 



lated on the basis of these points. Once the process moves from the geomet- 
rical data to point-based information, it can no longer regain the information 
about the geographical distribution of the location which was included in the 
25 parameter sets received from the mobile network. 

[07] The objective of the invention is to eliminate the drawback described 
above and to bring about a solution which enables the accurate location of the 
mobile through efficient usage of the information obtained from the mobile 
network. 
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Summary of the Invention 
[08] An objective of the present invention is to find a solution to improve the 

accuracy of the current location nnethods. Furthenmore, the objective is to 
achieve a solution not requiring extensive prior labor for achieving the im- 
proved accuracy, such as field measurements or data collection. 

[09] The invention utilizes location-dependent parameters available from a 
mobile network for determining the location of the mobile terminal. As men- 
tioned above, the parameters relating to a single measurement are in this 
context denoted as a parameter set. These parameters indicate the geo- 
graphical areas within which the mobile most probably is located, i.e. how the 
location of the mobile is geographically distributed over the area of the sys- 
tem. 

[10] In order to improve the accuracy of the system, the invention utilizes the 
geographical information provided by the parameter sets by forming matrices 
which indicate the probability distribution of the location of the mobile. A ma- 
trix is formed for a parameter set received from the mobile network. Each 
element of the matrix corresponds to a certain geographical area and contains 
a value which indicates the probability that the respective mobile is within the 
said area. At least one matrix formed for a mobile is stored as history data to 
be used in connection with a subsequent parameter set received for the mo- 
bile. When such a parameter set is received, the matrix corresponding to the 
current parameter set is fomned or retrieved, and the elements of the matrix 
stored earlier are updated. In the updating process the effect of the movement 
of the mobile on the matrix elements is taken into account, the movement be- 
ing the estimated movement occurring between the measuring events corre- 
sponding to the respective parameter sets. The location estimate is then de- 
termined on the basis of the element values of the matrix corresponding to the 
current parameter set and the element values of the matrix with the updated 
values. 



[11] Thus the invention utilizes history data in the form of at least one previ- 
ously fonned matrix, which is updated according to the estimated movement 
of the mobile. By employing the matrices, the original data received from the 
mobile network on the geographical distribution of the location of the mobile 
can be retained. The accuracy of the system can hereby be improved by 
processing the current matrix by one or more previous matrices associated 
with the mobile in question. 

[121 In one preferred embodiment of the invention, map information is em- 
ployed in determining the element values of the current matrix. The map in- 
formation preferably includes infonmation about the type of surface in the geo- 
graphical area related to an element. The description can be simple, such as 
road, forest, river, etc. The element values of the current matrix are then 
weighted according to the type of surface in question. 

[13] In another preferred embodiment map information is employed in updat- 
ing the element values of the matrix formed in connection with a previous pa- 
rameter set. 

[14] In addition to the improved accuracy of the system, a further advantage 
of the invention is that it does not require any special mechanism or laborious 
steps for forming the matrices. The matrices can be fomried "on-the-fly" based 
on the parameter set received, or they can be formed in advance and stored 
in a database. In the latter case the parameter set received Is used as a 
search key for retrieving the correct matrix from the database. 

[15] A still further advantage of the invention is that it is not dependent on the 
network implementation but can be applied to any network where at least one 
parameter dependent on the location of the mobile is available enabling a ma- 
trix to be formed. The method can therefore be used on top of network imple- 
mentations based on different location-dependent parameters. 



Brief Description of the Drawings 



[16] In the following, the invention and its preferred embodiments are de- 
scribed more closely with reference to the examples shown in Figures 1 to 9 
in the accompanying drawings, wherein: 

illustrates the location-dependent information available from a typical 
cellular network utilizing omni cells, 

illustrates the location-dependent information available from a typical 
cellular network utilizing sectored cells, 

illustrates a preferred embodiment of a system in accordance with 
the present invention, 

illustrates a sample matrix formed on the basis of a parameter set 
received from the mobile network, 

is an example of a flow diagram illustrating the method of the inven- 
tion, 

is an example flow diagram illustrating the processing of a preceding 
matrix, 

illustrates how the matrix of FIG. 4 is changed when it is updated 
once, assuming that all directions of mobile movement are equally 
probable, 

illustrates how the matrix of FIG. 4 is changed when it is updated 
once, assuming that the mobile is traveling directly eastwards, and 
illustrates one way of forming the matrix which corresponds to a 
parameter set received. 

Detailed Description of the Invention 
25 [17] As mentioned above, the method of the invention applies to various 

kinds of location-dependent information. Depending on the particular cellular 
system, the location-dependent information provided by the network can be 
signal strength or signal delay, for example. The determination of the location 
in current cellular networks is widely based on the Timing Advance value, be- 
30 cause the Timing Advance value is directly available from the network. There- 
fore, Timing Advance is in this context used as an example of the location- 
dependent signal information available from the mobile network for location 
determination. 
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[18] As is known, Timing Advance indicates how far the mobile most proba- 
bly is from the base station. FIG. 1 illustrates the location dependent infomna- 
tion provided by a network with omni-directional base station antennas, 
whereas FIG. 2 illustrates the same in connection with sectored cell sites. The 
network typically provides the Timing Advance infomnation as the minimum 
and maximum distance from the antenna (Rminand Rmax), in which case the 
mobile temiinal is with a certain probability between these limits, i.e. the 
hatched area A in the figures forms the Timing Advance zone defined by said 
limits. In addition to the Timing Advance information, the network provides the 
cell identifier CID, which identifies the cell in which the mobile terminal is lo- 
cated. This information can be given as the coordinates of the cell site. The 
network further provides an identifier for identifying the mobile in question 
from among the other mobiles and a time stamp indicating the moment of lo- 
cation measurement. In the case of a sectored cell, the network also provides 
the sector infomiation. Thus for each location determination the network pro- 
vides a parameter set which commonly includes the following information: cell 
identifying data, such as the coordinates of the Base Transceiver Station, 
Timing Advance information, such as Rmax and Rmm, an identifier for identify- 
ing the mobile in question from among the other mobiles, a time stamp indi- 
cating the moment of the location measurement, and possibly the sector in- 
formation. 

[19] FIG. 3 illustrates the key elements of the system according to the pre- 
sent Invention. It is assumed here that the mobile network is a GSM Public 
Land Mobile Network. Communication between the network and a mobile termi- 
nal MS in a cell takes place via a radio path by way of a Base Transceiver Sta- 
tion (BTS) 31. The Base Transceiver Stations are connected to Base Station 
Controllers (BSC) 32. Several Base Transceiver Stations are usually under the 
control of one BSC, and several Base Station Controllers are connected to one 
Mobile Switching Centre (MSC) 33, which canies out the main switching func- 
tions of the mobile network. In addition, the MSC connects the mobile network 
with external networks. For positioning purposes, the MSC is connected to a 



Gateway Mobile Location Center (GMLC) 34, which collects mobile position- 
ing information in a positioning database 35. As mentioned above, it is as- 
sumed in this context that the GMLC supports Cell ID and Timing Advance. 
This means that for each location determination the GMLC provides a pa- 
rameter set, including the parameters discussed above in connection with 
FIG. 1 and 2. 

[20] The parameter sets available from the mobile network are processed in 
an accuracy server 38, which receives location requests from external ob- 
jects, such as service applications residing in the network to which the accu- 
racy server is connected. In order to be able to determine a location estimate 
on the basis of a parameter set received, the accuracy server forms a matrix 
on the basis of the parameter set or retrieves a matrix corresponding to the 
parameter set from among a plurality of matrices formed in advance and 
stored in a matrix database 39a. The accuracy server further uses mobile- 
specific history data stored in a history database 39b. This data comprises 
matrices obtained in connection with the parameter sets received earlier for 
the mobile in question. The content of these matrices is discussed below. In 
determining the element values of the matrices, the accuracy server prefera- 
bly further employs map information stored in a map database 39c, although 
the use of map infomiation is not necessary. The same applies to the use of 
in-advance calculated matrices, i.e. the matrix database is not necessary, but 
the matrix corresponding to the current parameter set can be calculated in 
real-time. 

[21] The information required by the accuracy server for processing the pa- 
rameter sets can also be stored in a single database. The information is pref- 
erably stored in connection with the accuracy server, although it can be dis- 
tributed in the network to which the accuracy server is connected. 

[22] Upon receiving a parameter set, the accuracy server calculates a matrix 
corresponding to this parameter set. On the other hand, if the matrices have 
been calculated in advance, the accuracy server retrieves said matrix from the 



database 39a, and the parameter set is used as a search key for finding the 
correct matrix. The matrix contains nxn elements, each covering a certain 
geographical area and each including a value indicating a probability that the 
mobile is within the area covered by the element. FIG. 4 shows an example of 
a matrix comprising 10x10 elements. In practice, the number of elements is 
much greater, e.g. 40x40 elements. Each element corresponds to a coordi- 
nate pair Xi/Yi (i=0...9), which defines the geographical area covered by the 
element. In practice, the size of each element can be 25x25 m, for example, 
in which case a 40x40 matrix would cover an area of one square kilometer. In 
the example of Fig. 4, the element values are given as percents. The ele- 
ments covered by the Timing Advance zone are given the value of 100, while 
the other elements are given the value of zero. 

[23] FIG. 5 illustrates the estimation process performed by the accuracy 
server. As mentioned above, upon receiving a parameter set the accuracy 
server calculates a matrix corresponding to this parameter set or retrieves 
said matrix from the database 39a using the parameter set as a search key 
(step 51). The parameter sets can go through a normalization process so that 
the parameter sets which equal one another with a given accuracy are 
mapped to the same matrix. If the accuracy server retrieves the matrix, the 
matrices are preferably formed using a suitable calibration mechanism, such 
as the GPS system. In other words, when the matrices are formed before the 
commissioning of the system, the distribution of the location of the mobile is 
determined by the calibration mechanism and the matrix is formed on the ba- 
sis of the results given by the calibration mechanism. A combinatory mecha- 
nism can also be used, for example, so that the matrices calculated in ad- 
vance are further processed In the accuracy server in order to obtain the final 
matrix describing the probability distribution. In a preferred embodiment of the 
invention, this is done by using map information to adjust the element values 
of the matrix. 
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[24] In this embodiment the accuracy server retrieves map information in 
order to weight the element values by this data. The matrices to be weighted 
can originate from a database or can be formed in real-time by the accuracy 
server. Map infonnation refers here to any infonnation which correlates with 
5 the location probability or with the movement probability of the mobile. The 
map information used by the accuracy server is preferably in the form of map 
matrices similar to the matrices formed on the basis of the parameter sets, 
except that in the map matrices an individual element indicates the type of the 
surface covered by the element. The classification into surface types can be 

10 simple, such as road, forest, river, etc., each type being assigned a weight 
value. Some types are stronger than others; for instance, if a road goes 
through a forest, the subscriber is more likely traveling along the road rather 
than through the forest, so the type of the element is set to "road". The map 
matrix can be formed in real-time from the information retrieved from the map 

15 database, or the infomnation in the database can be in the form of matrices. 



■^i [25] The element values of the matrix formed on the basis of the parameter 

set received are then weighted according to the type of land in question (step 
53), whereby a weighted matrix is obtained. A previous matrix obtained in 

^ connection with a preceding parameter set relating to the same mobile is also 

20 retrieved from the history database 39b, and this matrix is supplied to an up- 

il date process. In this context, the previous matrix is also called the history ma- 

trix. In the update process (step 54) the effect of the estimated movement of 
the mobile on the location distribution is taken into account by spreading the 
probability distribution of the previous matrix according to the movement data. 
25 This operation is based on the information available about the movement of 
the mobile, such as speed and/or direction information. Map correlation cou- 
pled with movement information can be used to yield more accurate results. If 
no information on the movement is available, the probability distribution is 
spread evenly in all directions. As a result of the update process, a spread 
30 history matrix is obtained. 



[26] The current matrix and the spread history matrix are supplied to a com- 
bination process where a combined matrix of the same size is determined ac- 
cording to predetennlned rules (step 55). At least one location estimate is 
then calculated on the basis of the combined matrix (step 56). Furthermore, 
the combined matrix is supplied to the history database, where It replaces the 
previous matrix, I.e. the next time the previous matrix is retrieved, the com- 
bined matrix calculated in connection with the then preceding parameter set is 
obtained as the previous matrix to be updated. 

[27] FIG. 6 Illustrates the updating of the previous matrix, in which the values 
of the previous matrix are updated according to the assumed movement of 
the mobile. The update process Includes one or more calculation rounds. Dur- 
ing each round the probability distribution earned by the history matrix is 
spread once. The process employs three parameters which define how many 
rounds are needed each time; n, T, and T1 , where n indicates the number of 
times the spreading has been performed, T indicates the time elapsed from 
the previous update process (i.e. the time elapsed since the location of the 
mobile was estimated previously), and T1 is a fixed time interval which is de- 
ducted from T every time the matrix is spread. 

[28] The process first forms n1 xn1 basic coefficients wij used in the spread- 
ing process (step 60). In this example n1=3 and all the coefficients get the 
value of 1 as follows: 
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[29] Parameter n is then given the value of zero (step 61 ), since no spread- 
ing has yet taken place, and the value of T1 is deducted from the value of T 
(step 62). 

[30] The basic coefficients are then adjusted according to the movement 
information relating to the mobile, such as the direction and speed information 
(step 63). Map information can also be utilized in detennining the adjusted 



coefficients. It is assumed in this example that the mobile is moving directly 
eastwards, whereby the adjusted coefficients v/-,,] are as follows: 
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[31] In a spreading process the adjusted coefficients are then used to adjust 
the value of each element ek,i of the history matrix as follows (step 64): 

1+1 7+1 

M 

where M is the sum of adjusted coefficients deviating from zero, and 
.Ekj is the element value of the spread matrix. Thus the spread matrix is ob- 
tained so that each element of the previous matrix is replaced by a weighted 
average of itself and those of its neighbors which do not get a zero weight. In 
practice the number of elements affecting the adjusted value of each element 
can be greater, for example 16 (4x4 weights) or 25 (5x5 weights). 
[32] After the first spreading round, the accuracy sen/er tests whether T is 

still greater than zero and n smaller than a predefined limit nmax (step 65). If 
this is the case, the above-described spreading steps 63 and 64 are repeated. 
At the beginning of each new round, the value of n is incremented by one 
(step 66), and the value of T1 is deducted from the current value of T (step 
62). A new matrix is calculated as long as T remains greater than zero, and n 
smaller than a predefined limit n 

max- 

[33] As to step 63, it is further to be noted that the adjustment of the basic 
coefficients by means of the movement information may be performed only 
once for each update process (i.e once for each step 54). The adjustment ac- 
cording to the map information has to be made separately for each matrix 
element during each spreading round, since the map information used to ad- 
just the coefficients is element-specific, i.e. depends on the area defined by 
the element. 



[34] FIG. 7 illustrates how the matrix of FIG. 4 has changed after one spread- 
ing round when the above basic coefficients are used. Thus in this case no 
direction-based weighting is employed, but it is assumed that all directions are 
equally probable. As can be seen from the figure, the probability distribution 
has spread; the number of elements having a non-zero value has increased 
from 17 to 45. 

[35] FIG. 8 illustrates how the matrix of FIG. 4 has changed after one spread- 
ing round when the above-described adjusted coefficients are used. Thus, in 
this case it is assumed that the mobile is traveling eastwards. Now the num- 
ber of elements having a non-zero value has increased to 37. Utilizing the 
map information in the spreading process means that the probabilities spread 
faster along roads than along forests, for example. 

[36] The spread history matrix and the current matrix can be combined in 
various ways. However, there are always two probability values, A and B, cor- 
responding to a certain area, one from the element con-esponding to the re- 
spective area in the current matrix and another from the element correspond- 
ing to the respective area in the spread matrix. In the combination process 
(step 55), a new element value is determined for the respective element of the 
combined matrix based on these two values. There are at least three possible 
ways to determine the new element value for each element of the combined 
matrix: 

1 . calculate the sum of the two values, i.e. A+B, 

2. calculate a weighted sum of the two values, i.e. rxA+(1-r)xB, 

3. select the value which is the highest one, i.e. max{A, B}. 

As mentioned above, the combined matrix formed in one of the 
above ways is stored as a new history matrix. 

[37] Based on the values of the combined matrix, the final estimate(s) can 
also be determined in various ways. At least the following ways are possible: 

1. The coordinates of the element containing the highest value are 
selected as the location estimate. 



2. The coordinates corresponding to the center of gravity of the ele- 
ment values are selected as the location estimate. The center of gravity here 
refers to an average (calculated in a desired way) of the coordinate values 
associated \mth the element values. 

3. A certain area, which is obtained based on the previous estimate 
and the speed and/or direction information, is examined. The element contain- 
ing the highest value within that area is then selected as the location estimate. 
Another alternative is to select as the location estimate the coordinates corre- 
sponding to the center of gravity of the element values of the area. 

4. Local averages can be calculated for element groups of the de- 
sired size. For example, averages can be computed over local areas of 
100x100 m^, and the center point of the local area having the highest average 
can be selected as the local estimate. The form of the local area can also be 
other than a square, for example a circle. 

5. Areas of predetermined form, such as circles, are drawn around 
estimate candidates, each area covering a certain proportion of the total sum of 
the elements. The candidate is then selected for which the said area is the 
smallest. 

[38] The best alternative can be selected by testing the system in the desired 
environment and finding the alternative yielding the best results. 

[39] The formation of the matrix corresponding to the current parameter set 
(step 51) is discussed next. As mentioned above, the matrix can be formed 
directly in real-time on the basis of the parameter set, or it can be retrieved 
from among the matrices produced in advance by using the parameter set as 
a search key. A third alternative is to use a combination of these two methods. 

[40] When the matrix is formed in real-time from the parameter set, the geo- 
metrical counterpart of the parameter set, i.e. the area defined by the parame- 
ter set, can be superimposed on an empty matrix covering the area in ques- 
tion. FIG. 9 illustrates this by showing a Timing Advance zone 90 superim- 
posed on a 7x7 matrix.. Elements are then assigned values which reflect cor- 
respondence between area represented by the element, and the area that is 
actually covered by the parameter set. Thus for example. Fig. 9 the shape 90 
represents the geometrical area indicated by the parameter set. In the case of 



element 93 only about 40% of the element area is covered by shape 90, and 
thus the element is assigned the value 4. while element 95 is covered by 
shape 90 to about 70% of its area, and thus is given the value 7. As the ma- 
trix values range from one (no coverage) to ten (full coverage), these values 
5 can then be adjusted employing the map information. Alternatively, all the 

elements overlapping the geometric counterpart can be assigned the value of 
one, while the rest of the elements, which do not overlap with the geometrical 
counterpart, are assigned the value of zero. 

[411 When the matrix is retrieved from a database containing matrices 
10 formed in advance, it is preferable to measure the actual probability distribu- 
tion by means of a suitable mechanism, such as the GPS. Thus, in this case 
the matrix elements indicate the probability distribution measured in advance 

O in real environment. The parameter sets which are substantially equal to each 

£3 

other, i.e. which with a desired accuracy have the same parameter values, are 
15 associated with the same matrix. 

-1 [42] The third alternative which combines the above mechanisms could be 

such that an in-advance calculated matrix is always used when it is available, 
f ^ while the matrix is computed in real-time only when no calculated matrix is 

^ available in advance. This is because a matrix determined in-advance proba- 

y 20 biy represents a more reliable perception of the probability distribution. 

[43] Although the invention was described above with reference to the exam- 
ples shown in the appended drawings, it is obvious that the invention is not 
limited to these, but may be modified by those skilled in the art without depart- 
ing from the scope and spirit of the invention. In a simplified form of the inven- 

25 tion, no map information or direction information is used, i.e. the probability 

distribution is spread evenly in all directions. More than one history matrix can 
be employed for processing the current matrix. The network elements used 
can also vary. For example, the parameter sets can be received from any 
network entity having access to them, and the steps of the method can be 

30 distributed among different network elements. In the future a mobile may per- 
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form some or all of the steps of the method. Depending on the content of the 
parameter set, a mobile terminal might simply need external information re- 
quired to form the matrices in order to be able to perform the steps of the 
method. In a further embodiment, the mobile terminal can download prede- 
termined matrices from the network, i.e. matrix formation is not necessarily 
needed in the terminal. Upon receiving a parameter set, the mobile terminal 
retrieves or downloads from the network an in-advance calculated matrix 
which corresponds to the parameter set received. In this respect the mobile 
terminal can thus act similarly as the accuracy server. 



